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We present predictions for the formation of (anti)nuclear 
bound states in nucleus-nucleus reactions at RHIC energies. 
The phase space coalescence method is used in combination 
with RQMD-v2.4 transport calculations to demonstrate the 
relevance of particle production as well as the longitudinal 
and transverse flow components. The formation of deuterons 
follows an approximate scaling law proportional to the rela- 
tive freeze-out densities of nucleons and produced secondaries. 
For antideuterons, an additional suppression appears that is 
proportional to the number of nucleons, pointing toward mul- 
tiple rescattering and absorption prior to freeze-out. 



Nuclear clusters have been a useful tool to establish 
collective effects throughout the history of heavy ion re- 
actions. Their production rates have provided evidence 
for low temperature phase transitions [Q, their spec- 
tral distribution shows particular sensitivities to collec- 
tive flow transverse expansion |^J^-^| and poten- 
tial forces @1(J. In case of strong enough "cooling" 
of the emitting source and collective motion even the 
study of bound states with a considerable fraction of an- 
timatter [ p"l| , strangeness or even charm |l^ | becomes 
possible. Light antimatter clusters up to A = 3 have 
already been identified |H5), while the search for states 
with strange constituents is ongoing. Deuterons and an- 
tideuterons are the simplest composite objects and are 
useful in establishing expansion and correlations in the 
emitting source. Volume-expansion due to secondary in- 
teractions tends to diminish the cluster yields as par- 
ticle production rises both with the beam energy and 
the system size [jl4|. Counterbalancing effects can be 
expected from collective flow components that increase 
cluster multiplicities and reduce effective source radii as 
compared to the actual size of the system 0. It should 
be emphasized that the predictions of different transport 
models for collisions at RHIC energies already offer large 
differences in the most basic observables. For example, 
the total number of pions predicted at midrapidity varies 
by factors «2 in comparisons with parton cascade and 
RQMD-type calculations Jig] . Constraints from nuclear 
bound state analyses (i.e. fragment production) should 
complement those from single inclusive hadron spectra 



and pion/proton interferometry in order to distinguish 
the different model scenarios. Moreover, collective mo- 
tion, temperatures and position densities are related to 
entropy production and pressures, major assets in the 
search for QCD-phase transitions. We therefore suggest 
cluster analyses will remain a relevant tool in the up- 
coming experiments at RHIC. In this paper we present 
predictions for (anti)proton and (anti)deuteron produc- 
tion based on the transport approach RQMD v2.4 |ilf | 
in combination with a phase space coalescence frame- 
work 010 • Thc RQMD-transport model is a semi- 
classical microscopic approach which combines classical 
propagation with stochastic interactions. Color strings 
and hadronic resonances can be excited in elementary 
collisions. Their fragmentation and decay lead to par- 
ticle production. Overlapping strings may form ropes, 
chromoelectric flux tubes with charges in higher dimen- 
sional representations of color SU(3). RQMD is a full 
transport theoretical approach to reactions between nu- 
clei and elementary hadrons. Since the model does not 
include light cluster production it is supplemented with a 
coalescence "afterburner." The coalescence is performed 
by projecting the 2-body phase space density given by 
the microscopic transport onto bound state wave func- 
tions in Wigner-space. This method has been success- 
ful in the description of deuteron formation and proton- 
deuteron correlations at lower beam energies |7|,p|,|l7| . 
The basic observables that demonstrate phase space sen- 
sitivities are the rapidity and transverse momentum dis- 
tributions. Figure 1 (upper panels) shows predictions 
for (anti)protons and (anti)deuterons in central Au+Au 
collisions at full RHIC energy^]. Fig. 1 (lower panels) 
presents the rapidity dependence of average transverse 
momenta. The values exhibit approximately a factor 
of two higher transverse momentum in the composite 
objects. The average transverse momenta of antipro- 
tons (deuterons) are slightly larger than those of their 
matter counterparts. The predictions for particles and 
antiparticles shown in Fig. 1 are summarised in Ta- 
bles 1 and 2. The presence of flow components can be 



1 100A ■ GeVAu+ 100 A -GeVAu, impact parameters < b < 
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demonstrated by comparing with calculations where the 
freeze-out correlations of positions and momenta have 
been removed by hand via randomization (crosses on 
Fig. 1). As a result of the randomization process the 
total number of clusters at midrapidity is suppressed by 
orders of magnitude. The average transverse momenta 
at midrapidity in the composite spectra drop by approx- 
imately 30%, consistent with simple momentum coales- 
cence results. The collective component in the transverse 
momentum of deuterons and antideuterons is correlated 
with the number of scatterings: Collective expansion is 
highest in the midrapidity region while it becomes less 
pronounced close to projectile and target rapidities. A 
more detailed inve stig ation of transverse flow issues can 
be found in Ref. fig] . In the following, ratios of par- 
ticle yields are examined as measures of collectivity and 
(relative) freeze-out densities in longitudinal momentum. 
The particle yields show relatively small variations in ra- 
pidity space and the system shows strong evidence of 
a "collective" evolution towards freeze-out. In this case, 
the nuclear phase space densities can be approximated by 
average position densities and local momentum fluctua- 
tions, usually addressed as "temperatures" . For a purely 
thermal source the Nd/ Ap-ratio reflects the nuclear den- 
sity Nd/Np cx Np/V ~ (p p ) (where Nd is number of 
deuterons, N p is number of protons, p p is proton freeze- 
out density, V- freeze-out volume) and should not depend 
on the particle velocities. The transport calculations 
strongly deviate from such a scenario (see Fig. 2, a and c) 
showing Nd/Np- and Nd/N 2 - ratios that strongly vary as 
function of rapidity. This behaviour can be traced back 
to strong longitudinal flow components that lead to a par- 
tial separation of sources in beam direction. The effective 
"volume", V p is proportional to the ratio, oc p 2 /d, but 
seems to be enlarged proportional to the number of sec- 
ondaries N sec (here N sec is a number of produced mesons: 
N sec = iV(7r ± '°) + A^A'^ )) thus supporting the role of 
particle production and rescattering. This can be demon- 
strated by regarding the "scaled" ratio (Nd/N 2 ) * N sec 
for which the differences basically vanish. In semicen- 
tral reactions less expansion is found in the Nd/N 2 ratio 
while the scaled ratio turns out to be rather similar. The 
comparison of reactions with different projectile or target 
size has been found useful to assess the relative strength 
of hadronic expansion at SPS-energies flij| . Such size de- 
pendences are not seen in our centrality dependence anal- 
ysis although the effective nuclear volume seems to differ 
by factors 3-5. Note that the naive density-interpretation 
of the Nd/Np ratio can be somewhat flawed by trans- 
verse and directed flow correlations which change with 
rapidity. Antideuterons and antiprotons reveal differing 
emissions patterns that are dominated by the source cen- 
ter. The Nd/Np ratio (see Fig. 2, b) is rather flat and 
could lead to misleading conclusions since flow correla- 
tions appear to be close to none. The presence of flow ex- 
presses itself in the scaled and unsealed ratios N^/Np and 
(N^/Np) *N sec . They are only consistent with a volume- 



type scaling close to midrapidity but otherwise deviate 
considerably. This deviation can be explained by nuclear 
absorption which is larger due to the presence of higher 
nuclear densities and lower numbers of secondaries, par- 
ticularly in the domain \y\ > 3. As a consequence, the 
differences become most prominent in semi-central reac- 
tions with less scattering of baryons towards midrapidity. 
One further suggestion to address antimatter-absorption 
has been suppression at low transverse momenta as well 
as correlations in longitudinal and transverse flow (event 
plane asymmetries) [ pT| . In agreement with such a sce- 
nario, the average transverse momenta for antinucleons 
are slightly larger than those for nucleons. It has been 
suggested [|9| that the "homogeneity volume" deduced 
from cluster production should be compared to radius 
parameters extracted by HBT-type two-particle correla- 
tion analysis. We would like to point out that such an 
analysis suffers from major uncertainties which are re- 
lated to the widths of nuclear clusters in position and 
momentum space. Unlike HBT-type analysis, where the 
correlation strength can be scanned as a function of the 
relative momentum of particle pairs (qi,q s ,qo)i the pro- 
jection on boundstate wave functions involves both in- 
tegrations over relative momenta and positions. Hence, 
cluster analyses are more suitable to study the average 
phase space volume. Details of the emitting source in 
position space cannot uniquely be addressed unless the 
size and shape of the (local) momentum fluctuations are 
known. A solution to this caveat could be the study 
of various cluster types such as deuterons and He with 
wave functions (correlations) very different in position 
and momentum space. Heavier clusters, in addition, pro- 
vide less sensitivity to local momentum fluctuations than 
the loosely bound deuteron state. As collective flow be- 
comes strong enough, the characteristic scaling of the 
relative yields can give access to the flow and density ge- 
ometry ("event shape") [po[ . Further insight can be ex- 
pected from proton-proton and proton-deuteron correla- 
tion analysis (see for example |17f| ) which should confirm 
the shape of deuteron phase space densities. In summary, 
using the transport model RQMD(v2.4) and a coales- 
cence afterburner that projects RQMD's two body phase 
space densities onto fragment wave functions in Wigner 
space, we study the production of nucleons and deuterons 
(and their antiparticles) for the central 100A-GeV Au 
+ 100A-GeV Au collisions. Rapidity distributions and 
average transverse momentum exhibit strong longitudi- 
nal and transverse flow components. As a consequence, 
composite ratios are closely related to the position space 
distributions of the nucleons and produced hadrons close 
to freeze-out. Deuteron formation is consistent with a 
scaling relationship proportional to the relative densities 
of nucleons and secondaries at similar rapidities. The 
spectra of antideuterons are strongly modified during the 
course of the reaction. Contributions from antimatter 
absorption lead to deviations of antidcutcron production 
from the naive volume scaling and slightly higher trans- 
verse momenta as compared to their matter counterparts. 
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Yet to be addressed are observables such as elliptic [ pi|j22| 
or directed [0j2^| flow patterns which are more sensitive 
to details in the early and late event shape. We are grate- 
ful for many enlightening discussions with Drs. S. John- 
son, D. Keane, S. Pratt, H.G. Ritter, S. Voloshin. This 
research used resources of the National Energy Research 
Scientific Computing Center. This work has been sup- 
ported by the U.S. Department of Energy under Contract 
No. DE-AC03-76SF00098 and W-7405-ENG-36, the En- 
ergy Research Undergraduate Laboratory Fellowship and 
National Science Foundation and the Marie-Curie Re- 
search Training Grant No. FMBICT961721. 

t B. Monrcal is at Lawrence Berkeley National Labo- 
ratory through the Center for Science and Engineering 
Education. 



[2. 

[3] 

[4] 

[5, 
[6 

ir 

[8] 

[9] 

[10 
[11 

[12 



[13. 
[14 

[is; 



[16. 
[17 



nucl-ex/9805006 



subm. 



J. P. Bondorf, R. Donangelo, I.N. Mishustin, H. Schultz, 
Nucl. Phys. A444, 460 (1985); J. Pochodzalla et al, 
Phys. Rev. Lett. 75, 1040 (1995) 

M.A. Lisa et al, (EOS-collab.), Phys. Rev. Lett. 75, 2662 
(1995); M.D. Partlan et al, (EOS-collab.), Phys. Rev. 
Lett. 75, 2100 (1995) 

L. Ahle et al, (E802-collab.), Phys. Rev. C 57, 1416 
(1998) 

J. Barrette et al, (E877-collab. 
to Phys. Rev. C 

K.S. Lee, U. Heinz, Z. Phys. C48, 525 (1990) 

I.G. Bearden et al, (NA44 Collab.), Phys. Rev. Lett. 78, 

2080 (1997) 

R. Mattiello, A. Jahns, H. Sorge, H. Stoecker, W. 
Greiner, Phys. Rev. Lett. 74, 2180 (1995); R. Mattiello, 
H. Sorge, H. Stoecker, W. Greiner, Phys. Rev. C55, 1443 
(1997) 

A. Polleri, J.P. Bondorf, I.N. Mishustin, Phys. Lett. 
B419, 19 (1998) 

J.L. Nagle, S. Kumar, D. Kusnezov, H. Sorge, R. Mat- 
tiello, Phys. Rev. C53, 367 (1996) 

P. Danielewicz and Q. Pan, Phys. Rev. C46, 2002 (1992) 
M.Bleicher, C. Spieles, A. Jahns, R. Mattiello, H. Sorge, 
H. Stoecker, W. Greiner, Phys. Lett. B361, 10 (1995) 
J. Schaffher-Bielich, C. Greiner, H. Stoecker, A. P. Vis- 
cher J. Phys. G23, 2107-2115 (1997) ; J. Schaffner- 
Bielich, A.P. Vischer, Phys.Rev. D57, 4142 (1998) 
Quark Matter '96 proceedings, Nucl. Phys. A610, (1996) 
H. Sorge, J.L. Nagle, B.S. Kumar, Phys. Lett. B355, 27 
(1995) 

Talks presented on the Workshop on Transport Models at 
the RHIC98 Summer School, Brookhaven National Lab- 
oratory 

H. Sorge, Phys. Rev. C52, 3291 (1995) 
S. Panitkin, (E877 Coll.) , Proc. Int. Conf. Heavy Ion 
Physics at the AGS 96, eds GA. Pruneau et al, WSU- 
NP-96-16, 147 (1996) 



[18] B.Monreal, et al. |aucl-th/9990408| 3. 

[19] W.J. Llope et al., Phys. Rev. C 52, 2004 (1995); R. 

Scheibl and U. Heinz, PRC 59, 1585 (1999). 
[20] A. P olleri, R. Mattiell o, I.N. Mishustin, J.P. Bondorf, e- 

print |nucl-th/9904012| , subm. to Phys. Lett. B 
[21] J.Y. Ollitrault, Phys. Rev. D48, 229 (1992); Phys. Rev. 

D48, 1132 (1993) 
[22] H. Sorge, Phys. Rev. Lett. 78, 2309 (1997) 
[23] H. Appelshaeuser et al, (NA49-collab.), Phys. Rev. Lett. 

80, 4136 (1998) 



Rapidity 


P 


d 


P 


d 


0.25 


15.6 


0.042 


7.4 


0.0092 


0.75 


15.9 


0.049 


7.1 


0.0098 


1.25 


16.2 


0.057 


7.0 


0.0109 


1.75 


16.9 


0.078 


6.6 


0.0093 


2.25 


18.9 


0.098 


5.9 


0.0095 


2.75 


22.0 


0.150 


5.1 


0.0068 


3.25 


22.9 


0.195 


4.5 


0.0064 


3.75 


21.8 


0.202 


4.1 


0.0057 


4.25 


17.5 


0.208 


2.9 


0.0049 


4.75 


11.3 


0.159 


1.0 


0.0012 



TABLE I. Predicted production yield for protons, 
deuterons, antiprotons and antideuterons as a function of ra- 
pidity. 
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TABLE II. Predicted mean transverse momenta (in 
GeV/c) for protons, deuterons, antiprotons and antideuterons 
as a function of rapidity. 
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FIG. 1. Predicted rapidity distributions of (a) protons 
(filled circles), deuterons (open circles) and (b) anti-protons 
(filled squares) and anti-deuterons (open squares) for 
Au(200AGeV)Au reactions (b<3fm). (c) and (d) show 
the corresponding mean transverse momenta per nucleon 
as function of rapidity for protons, deuterons, antiprotons 
and antideuterons. Crosses indicate deuterons from calcula- 
tions without position-momentum correlations in the nuclear 
source. 
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FIG. 2. Predicted ratio of particle yields as function of ra- 
pidity for central (solid symbols) and semicentral (b=6fm) 
events (open symbols): (a)-(b) Nd/N p and N^/Np ratios, 
(c)-(d) scaled ratios {N d /Np)*N S ec and (N s /Np)*N sec and 
(e)-(f) iVd/iVp and N 3 /N? ratios. 
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